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The dynamics of a single hole (or electron) in the two dimensional Hubbard model is investigated. 
The antiferromagnetic background is described by a Neel state, and the hopping of the carrier is 
analyzed within a configuration interaction approach. Results are in agreement with other methods 
and with experimental data when available. All data are compatible with the opening of a mean 
field gap in a Fermi liquid of spin polarons, the so called Slater type of transition. In particular, this 
hypothesis explains the unusual dispersion relation of the quasiparticle bands near the transition. 
Recent photoemission data for Ca2Cu02Cl2 are analyzed within this context. 
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The nature of the insulating phase of the Hubbard 
model in a bipartite (square) lattice has attracted a great 
deal of attention. A clarification of the universality class 
of the transition will greatly help in understanding the 
nature of the "metallic " , gapless, phase away from half 
filling. A detailed discussion of the issue can be found 
in jy . Opposing viewpoints on the problem are presented 
in HH . An appealing approach is the scaling analysis 
discussed in [El. It is assumed that physical quantities 
can be expressed in terms of a single energy scale, A, 
which vanishes at the transition. A wealth of numerical 
results suggest that such scaling does exist. The data 
implies that typical length scales, £, scale as Ai. This 
implies the existence of a quasiparticle band, in the gap- 
less side, with a dispersion relation e/c cx k 4 |||| . 

The understanding of the Mott transition as function 
of filling in the Hubbard model in two dimensions (2D) 
requires a knwoledge of the dynamics of carriers intro- 
duced into the undoped antiferromagnetic (AF) system. 
This problem has been studied in the related t—J model 
HQ . The bands obtained by numerical methods are well 
described within the Born approximation ||], which ac- 
counts for hopping in one sublattice only. It has been 
shown that higher order corrections are small 111. The 
bandwidth scales, approximately, with J = 4t 2 /U. This 
can be understood from the fact that a spin flip is re- 
quired to restore the AF background to its original tex- 
ture after a hole hops between equivalent sites. Numer- 
ical results for the Hubbard model Jll 12\, show similar 
features, with a narrow band which is very flat between 
the (0, 7r) and (ir/2,ir,2) points of the Brillouin Zone. 
This band is also reasonably described by a generaliza- 
tion of the Born approximation to the Hubbard model 
p3| . As for the t—J model, this picture suggest that 
numerical results are well described in terms of dressed 
holes which hop within a given sublattice. 

The same picture is obtained from mean field calcula- 
tions. It has been shown that single holes in the Hubbard 
model induce the formation of inhomogencous spin tex- 



tures [Q, which can be described as spin polarons, or 
spin bags fl5| . An improved wavefunction results from 
hybridization of polarons localized at different sites, so 
that the translational symmetry is restored (ll| . This is a 
standard procedure in molecular physics, where Hartree- 
Fock (HF) solutions are improved by the Configuration 
Interaction (CI) method. The resulting polaron band is 
in reasonable agreement with other calculations Ijjjjl^] . 

In the following, we analyze the implications of the spin 
polaron band of the simple Hubbard model for the nature 
of the Mott transition. We also extend these calculations 
to include second and third nearest neighbor hoppings, 
in order to compare directly to experiments (l^^^]. Our 
results suggest that a mean field picture is a reasonable 
starting point for the analysis of the dynamics of holes in 
the Hubbard model in 2D, providing an intuitive picture 
of the Mott transition, and reproducing adequately the 
main features observed in the experiments. 

The usual HF approximation to the Hubbard model 
with nearest neighbor couplings only in a square lattice, 
gives an AF ground state at half filling. This solution re- 
produces the existence of a gap in the charge excitations, 
which scales with U , for large U. The approximation 
misses the low energy spin waves, which can be incorpo- 
rated by including the low amplitude fluctuations around 
the local minimum defined by the HF solution, that is, 
the Random Phase Approximation (RPA). The RPA re- 
produces the correct result in the large U limit |]l5| , ^l| . 
Thus, the mean field solution is a good starting point for 
the analysis of the hole dynamics. 

A naive generalization of the AF solution to the doped 
case can be obtained by filling the lowest lying states in 
the upper Hubbard band. Fig. [j] shows this (HF) band 
in the extendend Brillouin Zone, for t = 0.35 eV and 
U = 3.92 eV, in order to compare to experiments (see 
below). Due to the AF background, the bands have an 
additional periodicity, so that the point (%, ir) is degen- 
erate with (0,0). As the Hubbard model in the absence 
of other hoppings has electron-hole symmetry, the same 
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band describes electrons. The HF band shows a flat dis- 
persion along the (0, n) — > (tt/2,w/2) direction. The HF 
band can be written in the extended zone scheme as: 

A [A 2 

e k = — - y— + 4t 2 [cos(k x ) +cos(k y )} 2 (1) 

where A = y(n-f — n i)- Near the upper edge, this 
band can be approximated by: e k w (4t 2 /A)[cos(fe x ) + 
cos(fcj,)] 2 . It describes a hole with effective hoppings 
t 2 / A to the neighbor in position (2,0), and 2t 2 /A to 
that at (1,1). Hopping between different sublattices is 
not allowed. Around the point (Q, 7r), we can expand: 
e k w {t 2 / A)(k 2 — k 2 ) 2 . This quartic dispersion explains 
well the findings reported in H|| . 

The integrated density of states is dominated by the 
flatness of the dispersion around (0, 7r). Integrating 
over k, we find, in the Hartree Fock case, D(u>) ~ 
\og[t 2 1 (Uu>)]/ \ftjj near the upper band edge, and a similar 
dependence for the polaron band. 

The state c^l^o), where |<Fo) is the Slater deter- 
minant which describes the AF insulator, has negative 
compressibility and is unstable against charge inhomo- 
geneities JL4|. This topic was first discussed in [p2[ . As 
the compressibility, k, is proportional to (dfi/dn)^ 1 , it 
will be negative if the chemical potential is lowered upon 
doping The states near the gap edges are almost 

perfectly localized in one sublattice or the other. Thus, 
the additional electrons reduce the value of (n-f — ti]) by 
an amount directly proportional to the density. The gap 
is lowered by ^p, where n is the hole density per lattice 
site. Because of the singularity in the density of states 
at the gap edge, the difference between the chemical po- 
tential and the gap edge increases as {/i«n 2 , neglecting 
logarithmic corrections. This dependence cannot com- 
pensate the decrease in the gap, cx n, and the chemical 
potential is lowered upon doping. The transverse spin 
susceptibility also shows an instability as a consequence 
of the singularity in the density of states. Note, however, 
that the flatness of the band, which is a typical precursor 
of ferromagnetism J24j, does not imply a ferromagnetic 
instability. The interaction between carriers of opposite 
spin is greatly reduced, as they are localized in differ- 
ent sublattices. This analysis is in agreement with gen- 
eral studies of the possibility of phase separation in the 
Hubbard model p2[| . Similar results can be rigourously 
demonstrated in the limit of infinite dimensions p5|| . 

The previous shortcoming of the homogeneous mean 
field solution of the doped Hubbard model can be over- 
come by considering generic charge and spin textures [fl4[ . 
Among the many different solutions which are extrema 
of the Hartree Fock solutions, spin polarons and domain 
walls tend to be the most stable. We now analyze in de- 
tail the spin polaron, which is more stable at large values 
of U/t, and its energy is further lowered by derealiza- 
tion effects (see also below). The spin polaron, in the 



large U limit, can be approximately described as a self 
localized electron (or hole) for a cluster which contains 
at least five sites, and with a ferromagnetic alignment of 
the spins in its interior. The overlap between this so- 
lution and homogeneous solutions decays with the size 
of the cluster as 1/L 2 J26| , implying the formation of a 
localized state. This state is separated from the upper 
(lower) Hubbard band by an energy of order t in the 
large U limit. This splitting suppresses the contribution 
of the extra electrons to the spin susceptibility, curing 
the instabilities of the homogeneous HF approximation. 
In addition to this state, another localized level splits off 
from the lower (upper) band, and moves to an energy 
of order of a fraction of U from that band. Thus, the 
low energy spectral weight per electron (hole) is greater 
than one, in qualitative agreement with the arguments 
given in The Slater determinant which describes a 
solution of this type can be written as: 

1*0 = <T n * (E (e fo**^} i*°> ( 2 ) 

where |\&o) describes the antiferromagnetic insulator at 
half filling. Thus the wavefunction j^,-) is obtained from 
|\&o) by adding a given number of electron- hole excita- 
tions, plus an electron (hole) at site i. 




FIG. 1. Quasiparticle band structure for a single hole 
in the Hubbard band with t = 0.35eV and U = 3.92eV. 
Hartree-Fock approximation: dashed line. CI approximation 
(full line): dispersion relation, see Eq. (3), fitted to the nu- 
merical results for a 14 x 14 cluster. 

As discussed elsewhere [ pT|Jl8| , corrections to this inho- 
mogeneous HF solution can be of two types: low ampli- 
tude fluctuations, which can be studied within the RPA 
PH , and hybridization of solutions centered at different 
lattice sites |l6) . The combination of these solutions leads 
to wavefunctions of the type \^k) = e lkri and to 
the formation of a polaron band, whose width scales as 
t 2 /U in the limit of large U. 

The dispersion band for a 14 x 14 cluster and the same 
values of U and t used above, is shown in Fig. EJ. The 
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polaron (CI) band has the same shape than the HF band 
but is substantially narrower (around a 40%). We ad- 
scribe this narrowing to the fact that the polaron band is 
more strongly dressed by spin excitations (note, however, 
that the approximation used here only includes longitu- 
dinal, Ising like, modes). Because of the extension of the 
individual polarons, the band cannot be parametrized 
in terms of a few hopping parameters, as in the homo- 
geneous solutions described previously. The numerical 
results of Fig. [l| can be accurately fitted by, 

£k = £o + 4*n cos(fc x ) cos(ky) + 2t 2 o[cos(2k x ) + cos(2fcj,)] 
+4£ 2 2 cos(2fc a ) cos(2k y ) + * 3i [cos(3fec) cos(k y ) 
+ cos(k x ) cos(3k y )] + 2£4o[cos(4/c2,) + cos^fc^)]. (3) 

with t n = 0.130542 eV , * 20 = 0.06 2 56 eV, * 22 = 
-0.006130 eV, t 31 = -0.003963 eV, and * 40 = -0.000836 
eV. In the (0, it — rf) direction we can expand the disper- 
sion relation for small r\ as, 

£k ~ £0 + 4(-tn + <20 + *22 - 2t 31 + ) 

+277 2 (in - 2t 20 - 4* 22 + 10*3! ~ 8t 40 ) 
n 4 1 

+y (- gtii + 4*20 + 8*22 - 41< 3 i + 64* 40 ) . (4) 

Using the parameters given above we obtain ek f=a 
— 0.004r7 2 + 0.08?7 4 . This result shows that the quartic 
term dominates, as in HF and in agreement with |Q,|J. 
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FIG. 2. Spin (arrows) and charge distribution around a 
localized hole in a 14 x 14 cluster using the parameters from 
Ref. [20]. 

Experiments on the dynamics of holes in CuO planes 
suggest that a Hubbard model with nearest neighbor hop- 
pings only is insufficient |29|JT9| , |20| | . The inclusion of ad- 
ditional hoppings is straightforward within the scheme 
discussed here. We have calculated the electron and hole 
bands using the parameters given in [EG]: U = 3.92 eV, 
* = 0.35eV (J=0.125 eV), *' = -0.12cV and *" = 0.08eV. 
For hole doping this choice of parameters leads to the for- 
mation of well localized spin polarons at the mean field 
level. The spin and charge textures associated to the 
(hole) spin polaron are shown in pi, for a 14 x 14 cluster. 



We find no inhomogeneous solutions for electron doping 
(see below) 

The quasiparticle bands are shown in fig. j^] . The new 
hoppings break the electron-hole symmetry present in 
the simple Hubbard model and induce a quadratic dis- 
persion at the band edges, which dominates the quar- 
tic terms discussed so far. This dispersion relation im- 
plies that, upon doping, emptying the lowest lying states 
give a change in the chemical potential Sfi oc n that can 
be compensated by the decrease in the gap oc n. As 
a result, homogeneous mean field solutions can be sta- 
ble. The HF and CI bands have a similar shape, while, 
as in the case of the simple Hubbard model, the CI 
band is much narrower than the CI band (an approxi- 
mate reduction of 40% is noted). As remarked in ]2C| l 
the HF band is much wider than the one given by pho- 
toemission experiments. The CI band is in much closer 
agreement with the observed band (approximately 0.4 
eV wide) and with the numerical results for the corre- 
sponding t-J model reported in The shape of the 
hole band shown in Fig. 2 is also in full agreement with 
experiments. The CI band can be fitted with the disper- 
sion relation of Eq. (3). In this case the parameters are 
*n = -0.071428 eV, * 20 = 0.19 5 491 eV, * 22 = -0.016190 
eV, * 3 i = 0.003486 eV, and * 40 = -0.00 6 2 78 eV. The 
expansion near (0,7r) gives constant — 0.63?y 2 + 0.05?7 4 , 
which is mainly quadratic. 




(o,o) (o,tt) (Va.T/a) (0,0) 



FIG. 3. Quasiparticle band structure for a singly doped 
Hubbard model with second and third nearest neighbor hop- 
pings, using the parameters from Ref. [20] Hartree-Fock ap- 
proximation: dashed line (hole doping) and dotted line (elec- 
tron doping). CI approximation (full line): dispersion rela- 
tion, see Eq. (3), fitted to the numerical results for a 14 x 14 
cluster. 

A detailed analysis of the hole dispersion relation in the 
corresponding * — *' — *" — J model is presented in J30| . It 
is instructive to compare our results with those there. As 
mentioned at the beginning, the nature of the approxima- 
tion used here for the Hubbard model describes similar 
processes to those included in the Self Consistent Born 
Approximation used in p5| . However, when reducing the 
Hubbard hamiltonian to the *— J one, hopping terms oc J 
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are neglected. Thus, the bandwidth obtained from the 
t — J hamiltonian should be smaller than the exact. The 
approximation used here neglects the interaction of the 
hole with transverse spin waves, which probably leads 
to an overestimation of the exact bandwidth. Magnetic 
correlations around the hole are in qualitative agreement 
with (30) (sec fig.[@). 

For electron doping, instead, the homogeneous solu- 
tion is stable. The HF band is shown in Fig. The 
band has a significantly smaller effective mass near (0,7r) 
than in the case of hole doping, in agreement with [ fL9[ |, 
resulting in a larger quasiparticle peak as observed ex- 
perimentally [fj9| . 

The instability of the homogeneous mean field solution 
for the doped Hubbard model can also lead to the for- 
mation of stripes, or domain walls [[^[h] [33|. Similar 
solutions have been found in the t-J model [ |34p5| l). A 
transition from the undoped insulator to a doped sys- 
tem with stripes should be qualitatively different from 
the one described above. The stripes will remain, most 
likely, static, even after the inclusion of corrections be- 
yond HF. The lack of long range magnetic order can be 
understood from the weakening of the effective exchange 
coupling, and the enhancement of fluctuations |36|, as 
for the spin polaron solutions. Thus, we cannot rule out 
the possibility of a sharp discountinuity, as function of 
U/t and for a fixed low doping, between a gapless phase 
with delocalized spin polarons, and another with static 
domain walls. 

The above analysis shows that the properties of the 
insulating and the lightly doped phases of the Hubbard 
model can be reasonably understood within straightfor- 
ward extensions of mean field theory. The picture is sim- 
ilar to the spin bag model jD|[[7||8) . The main difference 
is that the spin bag approach uses homogeneous solutions 
as starting point. We have addressed the instabilities 
of these solutions, and showed how they can be over- 
come. The adequacy of mean field theory in describing 
the metal-insulator transition is due to the suppression 
of dynamic charge fluctuations by the formation of the 
staggered magnetization, as in infinite dimensions [j39| . 
The resulting picture leads to a flat quasiparticle band, 
in agreement with detailed numerical calculations. The 
density of states is singular at the band edges, enhanc- 
ing the effects of the interactions. A quartic dispersion 
relation is obtained around the point (0,7r). Including 
second and third nearest neighbors hoppings allows us 
to compare with recent photoemission data. This new 
hoppings enhance the quadratic terms in the dispersion 
relation and break the electron/hole symmetry. The CI 
hole band is much narrower than the HF band and in 
agreement (both shape and width) with the experimen- 
tal results. 
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